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Abstract. For the integrable iV-particle Calogero-Moser system with elliptic 



potential it is shown that the Lax operator found by Krichever possesses a classical 

r-matrix structure. The r-matrix is a natural generalisation of the matrix found 

OQ ■ 

recently by Avan and Talon ( |hep-th/9210128| ) for the trigonometric potential. The 

r-matrix depends on the spectral parameter and only half of the dynamical variables 



(particles' coordinates). It satisfies a generalized Yang-Baxter equation involving 
another dynamical matrix. 



1 Introduction 



• i 

The elliptic Calogero-Moser model is the system of N one-dimensional particles 
^ | interacting via two-particle potential 

V{q a p) = p(q a p), q a p = q a -qf3, (1) 

p being the Weierstrass function [J] with the periods 2uj\, luii- In terms of the 
canonical variables p a , q a (a — 1, . . . , N) 

{Pa,Pp} = {q a ,q/3} = 0, {p a , qp] = 5 a p (2) 
the Hamiltonian of the system is expressed as 

N 

w = Ei£ (3) 

(we deliberately omit a coupling constant which always can be removed by rescaling). 
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The Hamiltonian (|3|) with the potential (|I|) is known to be completely integrable 
0, |3|, H, [5|. The most effective way to construct the complete set of commuting 
integrals of motion is to represent them as the spectral invariants of a Lax operator 
- certain N x N matrix depending on the dynamical variables and, may be, on an 
additional parameter u called spectral. As shown in ||, for the commutativity of 
the spectral invariants tiL(u) n of the Lax operator it is necessary and sufficient that 
the Poisson bracket {L^(w), L^(v )} could be represented in the commutator form 

{L£(«Ua(«)} = E 

7172 

or, using the notation |7| 

L W = L^1, L (2) = 11®L, 



r^g 2 (u, v)L}\ (it) - m {u)r}\f 2 (u, v) 
-^{v,u)L%(v)+L^(v)r^(v,u) (4) 



as 

{L^(u),L^(v)} = [r^ 2 \L^\u)} - [r^ 21 \ L^{v)\ (5) 

where r^ 12 ^ is an N 2 x N 2 matrix depending, generally speaking, on the dynamical 
variables, and r^ 21 ) is 

r {21 \u,v) = Pr {l2 \v,u)P, (6) 

P being the permutation: Px <8> y = y ® x. 

In contrast with the well-studied case of purely numerical r-matrices 0, no 
general theory of dynamical r-matrices exists at the moment, apart few concrete 
examples and observations || |, |9], [U], [U], [12], |13[. Still, the collection of examples is 
rather sparse, and any new example of dynamical r-matrix could possibly contribute 
to better understanding of their algebraic and geometric nature. 

In the dynamical r-matrix has been calculated for an L-operator without 
spectral parameter corresponding to Calogero-Moser model with the trigonometric 
potential V(q) = 1/ sin 2 q or 1/ sinh 2 g. In [11| an elegant interpretation of the r- 



matrix was found in terms of hamiltonian reduction. However, the L-operator used 
in 



nj [11]] lacks the spectral parameter. Though the absence of the spectral parame- 
ter is not important for the construction of commuting integrals of motion, the spec- 
tral parameter is indispensable for integrating equations of motion and constructing 
the action-angle variables. The proper way of introducing a spectral parameter into 
L-operator was found by Krichever [TJJ . 

In the present paper the dynamical r-matrix is found for the Krichever's L- 
operator corresponding to the elliptical potential ([j]). The r-matrix found previously 
T0|1 is shown to be a degenerate case of our r-matrix. We derive also an equa- 



ni 



tion satisfied by the dynamical r-matrix which generalises the classical Yang-Baxter 
equation satisfied by numerical r-matrices. Our equation has a slightly more general 
form that the equation for a dynamical r-matrix proposed recently in flj| . 
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2 L-operator and dynamical r-matrix 



Krichever's L-operator [14| is the N x N matrix given by the formula 



Lp(u) = p a 8 a p + iQ(u, q a p){l - S a p), 



(7) 



or, in terms of the basic matrices Ep 



{Ep)p, - Saa'Spp', 



by 



N 



(8) 



a=l 



where Q(u, q) is expressed in terms of Weierstrass a functions 

a(u + q) 



Q(u,q) 



a(u)a(q) 



(9) 



For what follows, it suffices to know that a(z) is the entire function defined by 
the product M 



a(z) = z I] ( 1 - 



exp 



z 1 / z x 2 

— + 



OJr, 



2 W 



where ui mn = 2muji + 2nu2, 2ui >2 being a pair of periods. The ( and p functions are 
introduced through the derivatives of a 



CO) 



cr U 



o (z 



P(z) = -C(z) 



(10) 



The translations of z by 2ui } 2 and the reflection z — > —z act on the functions as 
follows 

a(z + 2u>i) = -a{z) exp[2(;z + uijCfa)], 

CO + 2u t ) = ((z) + 2C(w,), p(z + 2u{) = p{z). 

<T {- z ) = -a(z), C(-z) = -((z), p(-z)=p(z). (11) 
In the vicinity of z = the Weierstrass functions have the expansions 

a(z) =z + 0(z 5 ) 1 ((z) = z- 1 + 0(z 3 ), p(z) = z- 2 + 0(z 2 ). 

The spectral invariants trL(u) n , n — 1, . . . , N are independant commuting quan- 
tities 0. For instance, the total momentum V = Y.a=iPa an d the Hamiltonian 7i 
(0) are obtained from 



trL(u) = V, trL(w) 2 — H- V(u) 



(12) 
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the last equality following from the identity 



Q{u, qa/3)Q{u, qpa) = V(u) - V(q a/3 ) (13) 
which follows, in turn, from 

<j(u — v)a(u + v) 



<j 2 (u)a 2 (v) 



p(u) - p(v), (14) 
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Proposition 1 For the L-operator ([Jj the identity (fj^) holds with the following ma- 
trices r*- 12 ) et r*- 21 ) 

N 

r^ (u, v)=aJ2 Kl + £ ^< + £ d aP {E^ + E%) (15) 



where 



P Q l a 2 — p«l <> p«2 



-C(u - «) - C(v), c a/3 = r a f a = -Q(u - v, q af3 ), 



l aj3 



t Tb = r pp = -\Q(v,q a p), 



(16) 



-(21) 



being defined by (|6[) 



Proof. Notice, first of all, that, due to (|i~6|), the tensor 



P1P2 



has nonzero 



components only for combinations of four indices [sy^] having no more than two 
different indices. This property of r simplifies greatly verification of the identity (^) 
since, in this case, it suffices to consider only the combinations [Jgyj 2 2 ] in (|j) having 
no more than three different indices. The obvious symmetry of (|T6D with respect to 
the permutations of indices allows to reduce the task to the verification of (§) for 
only 14 combinations 



ill] 



ii]) [12L [21]) [ii]> [22]' [12L [2?]' [22]' [13]' [32]' [23]' [33]' [3?]' 



Substituting (|T5p into 
equalities 



ill] 
ill] 
ill] 
ill] 

ill) 
ill] 
ill] 
ill] 



{L\(u) 
{L\(u) 
{L\{u, 
{L\{u) 
{Ll(u) 

{L{« 
{L\{u) 
{Ll(u) 



L\{v)} 
L\{v)} 
L{(v)} 
Ll(v)} 
Ll(v)} 
Ll(v)} 
Ll(v)} 
Ll(v)} 



|4j) results for the above combinations of indices in the 
0, 

— L\{u)c2\{u — v) — a(v, u)Ll(v), 
a(u,v)Ll(u) + L\(v)c 2 i{v - u), 
c 12 (u - v)L\(u) + Lj(v)a(v, u), 
—Ll(u)a(u, v) — c 2 i(f — u)Ll(v), 
0, 

c 12 (u - v)(L 2 2 (u) - L\{u)) + c 21 (v - u)(L 2 2 (v) - L\(v)), 
0, 
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[{I}: {L\(u),Ll(v)} = 
[32] 

[Ml 
ill] 
Hi) 



ill] 



{Ll(u),Ll(v)} = 0, 

{L\{ulL\(y)} = d 13 (v)Ll(u)-d 12 (u)Ll(v), 
{Ll(u),Ll(v)} = -Ll(u)d 23 (v) + Ll(v)d 13 (u), 
{L\{u),L\{v)} = d 21 (v)Ll(u) + c 12 (u-v)Ll(u) 

+L\{v)d 13 {u) + Ll(v)c 31 (v - u), 
{Ll(u),Ll(v)} = -L\{u)d 23 {v) - L\(u)c 32 {u - v) 
-d 12 (u)Ll(v) - c 21 (v - u)L\{v). 



«l«2l 



Using then the substitution (Q) one obtains the trivial identities = for 
[11L [12]) [22]' [13]' an d [32]- After the substitutions (|T6|) the same trivial ide±±u±uj ±0 
obtained for = [%], and [*§]. 

In the same manner, the case [21] is reduced to the equality 

= (pi - p 2 ) [Q(u - v, q 12 ) +Q(v-u, q 21 )\ 
which follows from the identity 

Q(-u,-q)=-Q(u,q) (17) 

resulting from @ and ([II]). 

\ct1ct2~ 

= Q{v,q 21 )Q(u,qi 3 ) + Q{u - v,q 12 )Q(u,q 23 ) + Q(v,q 23 )Q(v - u,q 31 ) 
= Q(u,q 12 )Q(v,q 23 ) + Q(u,q 13 )Q(u - v,q 32 ) + Q(v - u, q 21 )Q(v, q 13 ) 



Analogously, for = and [l 3 ] one obtains, respectively, the identities 



which, after substitution of ([]), are reduced to the standard three-term quartic 
identity for u-function |L5] 



a(x — y)cr(x + y)u(z — t)a(z + t) 
+ cr(y — z)a{y + z)a(x — t)cr(x + t) 

+ a(z - x)a(z + x)a(y -t)a(y + t) = 0. (18) 

Similarly, the remaining cases [^^ 2 2 ] = [12], [ 2 \], and [22] lead, respectively, 
to the identities 

[C(v + q 12 ) - ((qi2)]Q(v, qn) = Q(u - v, q 2 i)Q(u, qi 2 ) + [((v - u) + ((u)]Q(v, q 12 ), 
[C(u + q 12 ) - C(<li2)]Q(u, qn) = Q(v - u, q 2 i)Q(v, q 12 ) + [((u - v ) + ((v)]Q(u, q 12 ), 
[t(v + ?2i) - C(Q2i)]Q(v, 921 ) = Q(u - V, qi 2 )Q(u, q 21 ) + [((v - u) + C{u)\Q{v, q 21 ), 
[C(u + qn) - C(qi2)]Q(u, qn) = Q(v - u, q 2 i)Q{v, q 12 ) + [((u - v) + ((v)]Q(v, q 12 ), 

which are reduced to another standard identity |15[ 

n n n \ cr(22;)cr(x + y)a(x — y) , , 

C x + ^ _ C x _ z + C (y _ z _ C y + z = \\ y ,(\s 7 r 19 

cr(x + z)cr(a; — z)cr(?/ + -2)cr(y — 2) 

being a derivative form of (|l8l) . 
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3 Degenerate cases 

The degenerate cases discussed below correspond to the infinite values of one or 
both of the periods lo\^. The corresponding potential V(q) becomes respectively a 
trigonometric/hyperbolic or rational function of q. According to []5| we shall mark 
these cases as I (V(q) = q~ 2 ), II (V(q) = sinh _2 g), and III (V(q) = sin~ 2 g), the 
case IV corresponding to the original elliptic potential V(q) = p(q). 

In the case II the periods are oj\ = oo, aj 2 = the Weierstrass functions 

degenerate into M 



sinh z exp 
and Q(u,q), see (|9]), into 



coth z — 



sinh 2 z ^ 3 



_ sinh( M + q) 
sinh u sinh q 



uq 

y 



(20) 



(21) 



One can simplify the corresponding L operator ([S]) removing the factors e u i a P/ 3 
by means of the similarity transform 



5 a/3 e 



uq a /3 



L(u) — ► L'{u) = lf(M)L(tt)^ 1 (tt), 

Calculating the Poisson bracket {L'^(u), L ,<y2 \v )} 

{L'M,L'W} = {W^L^ l \W^)-\W^L^\W^)- 1 } 
= W^W^{L^ l \ L^}(W^)-\W^)- 1 

+W^{W^,L^}(W^)- 1 L^(W^)- 1 

+iy (1) {£ (1) , w^}{w^)- 1 L^\w^y l 

-W {l) W {2) L^\W (2) )-\L^ l \W {2) }{W {2) )- l {W w )- 1 



(22) 



(23) 



(we have omitted the arguments u and v which are easy to restore) and using the 
identities 



W /(D{ jL (l) )W /(2) }(H /(l) ) -l (W /(2) ) -l 



- U -E 
3 

V -E 
3 



where 



N 



act 



a=l 



one can show that L' satisfies again the relation @ with the modified r matrix 
which is also described by the formulas flTSP and (|T^) where one should replace the 
Weierstrass functions with 



viz 



sinh z, 



C(z) 



coth z, 



P 



sinh z 



(24) 



and, respectively, 

_. , s sinhfw + q) 

Qiu, q) — > — : : — = coth u + coth q. (25) 

smh u smh q 

Quite analogously, the case III, corresponding to the periods uj\ = n/2 and 
oji = ioo, is reduced to the substitutions 

a(z) — > sin z, ((z) — > cot z, p — > — 

sin z 

and 

. sin(w + q) 

Qiu, q) — > = cot u + cot q. 

sin u sin q 

The case I corresponds to the periods lo\ = oo and u 2 = ioo and to the substi- 
tutions 

1 1 

a(z) — ► z, C(z) — ► -, p — > — , (26) 

z z z 

and 

Q(u,q)-^^- = - + -. (27) 
uq u q 

In [JTIJ the r matrices were calculated for the I, II and III type potentials and the 
L operators containing no spectral parameters. For the case I the L operator used 
in [ 10 1 has the form @ with Q(u,q) replaced by q^ 1 . The corresponding r matrix 



has the structure (|TJ) with 



11 

0, c a/ 3 — , d a p — — (28) 



It is easy to see that the above result can be reproduced from ours (p26|) , (P7|) by 
taking the limits u, v — > oo. The cases II and III need, however, some tricks. 

It is sufficient to consider the case II since the case III is obtained by the obvious 
change of the hyperbolic functions into the trigonometric ones. For the case II the 
L operator used in [[U| has the form (§) with Q(u,q) replaced by 1/ sinhg. The 
corresponding r matrix has the structure (fl5|) with 

a = 0, Ca/3 = - cothg a/3 , d a(3 = -77—^ (29) 

z smn q a f3 

In order to reproduce this result let us start with the L operator and the r matrix 
given by (§), ([15]) and (g§, and take the limit 

u — > +00, Q(u, q) - 



sinh q 

Removing the factors e qa/3 from Lg by means of the similarity transform 

L — ► WLW~\ Wj? = 6 af3 e- qa 

one obtains the L operator used in Jl0[. The modified r matrix is calculated then 
as previously ( ^3l) with the only difference that one has to coordinate the limits of u 
and v. The simplest choice is u — v = m/2. It allows to avoid the pole (u — v)^ 1 and 
to preserve the relation (||). The resulting r matrix coincides with the one found in 
10 . 



7 



4 Generalized Yang-Baxter equation 

It is well known that for a purely numeric r matrix a sufficient condition for the 
Poisson bracket defined by (|5]) to satisfy the Jacobi identity is the classical Yang- 
Baxter equation J7|. 

[ r (12) )r (13)] + [r (12) ;r (23)] _ [r (13) jr (32)] = g ( g 0) 

(in this section we omit, for the sake of brevity, the spectral parameters u, v, w 
corresponding, respectively, to the spaces 1, 2, 3 in the tensor product ® <g> 
C N ). 

The natural question arises of a proper generalization of (|30| ) for the dynamical 
r matrices. Starting with the Jacobi identity for the triple L^ 2 \ 

L®}, L®} + {{L< 2 >, L^}, L«} + {{L®, L™}, L®} = (31) 

and using (|5|) to estimate the Poisson bracket of two L operators we obtain the 
equality [|, |] 

[i?( 123 \ L«] + [i?( 231 U (2) ] + [^ (312) , ^ {3) ] = 0, (32) 

where 

^(133) = r (123) _ {r (l 3)) L (2) } + {r (l 2)) L (3) }; 

r (i23) De j n g kh e left-hand-side of fl30l). 
Trying to satisfy (^) with the ansatz 

£(123) = [X (123) ; L (2)] _ [y(123) ) L (3)] (33) 

we obtain the relation between X and Y 

y(123) _ (312) (34) 

or, at length, 

%X(V,») = X%£g(w,u,v). (35) 

We are led thus to the following generalization of the Yang-Baxter equation for 
dynamical r matrices 

j- r (12) j7 .13)l + |- r (12) jr (23)j _ ^(13)^(32)1 
_{ r (13) )L (2) } + {r (12) ;L (3) } 
_ [X (123) )L (2) ] + [ X (312) )L (3)] = (36) 

involving a new object — the X matrix. 

Notice that our ansatz (^) for R^ 123 ^ is more general than that proposed in [ 13| : 



£(123) = [X (123) ;L (2)_ L (3)] ) 

^(123) ^ging assumed completely symmetric with respect to the permutations of 123. 
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Proposition 2 For the L operator (Qj and the r matrix ( jT^ j the identity ( |5^[ ) /ioWs 
■ura£/i i/ie following matrices X^ 123 ) and X^ 312 ) 



X^ 23 \u,v,w) = -iJ2Q( w ^^) 

OL + P 



^111 

u rpaaa , x rp/3f3a . x pa/Ja , x 771/300 
-T g-^/3/3/3 4-^/3/3 4 " C/ /3a/3 



(37) 



where 



^1/32/33 - ^1 ^ ^2 ® ^3 ' 



X< 312 ) feeing given by jfflj). 



components only for combinations of six indices [^ 2 2 /3 3 3 ] having no more than two 
different indices, it suffices to consider only the combinations [a 1 ^ 2 ,?, 3 ] having no 



Proof. The verification of ( ftSD is quite similar to that of (^) in the proof of 
the Proposition 1. Again, one notices that, since the tensor X^ 1 ^ 2 ^ 3 has nonzero 

r«i«2"3l L^r^r 
L/3i/3 2 /3 3 

" "" L/3i/3 2 /3 3 

more than three different indices. The obvious symmetry of (^) with respect to 
the permutations of indices allows to restrict the list of the combinations to be 
verified with the total of 1 + 31 + 90 = 122 items, which is, however, too much to be 
done by hands. Fortunately, a computer-assisted calculation shows that 102 of 122 
equalities result, after substitutions flD, (|l^) and fl37|), in the trivial identities = 0. 



It remains to verify only quite manageable amount of 20 identities corresponding to 
the combinations = [}»], [\\% [ft], [&], [%]], [}»], [%% [£»], [?£], [J 22 ], 

[121] [1121 [122] [1221 [1121 [1211 [1231 [1231 [1231 [1231 
1.122b Ll22j; L221J5 [212J5 Ll23j) 1. 132b Lll2j) Ll3lJ) L 23lJ ' L312J - 

These 20 identities can be divided into 3 groups. The first group includes 6 
identities corresponding to the combinations [HI], [\H], [\H], [\ll], [III]: [312] which 
are analogous to the identities [gf] and [33] (see the proof of the Proposition 1) and 
are reduced to the standard identity (111 



The second group includes 12 identities corresponding to the combinations [\l\], 

[111], [ft], [ift, [lill W, [ml> Will W, [ml> [22IL [SI, which, like the identities 
[ft, [21], [n], and [33] from the quoted proof, are reduced to the second standard 
identity (0). 

The remaining 2 combinations [ 1x2] j [121] lead, respectively, to the identities 

C(u - v ) + ((v) - C(u - w) - ((w) = 

Q(u - v, qu)Q(u - w, ggi) - Q(v, q 21 )Q(w, q 12 ) 
Q(v - w,q 2 i) 

C(u -v) + C(v) - C(u -w)- C(w) = 

Q{u - v, q 21 )Q{u - w, q 12 ) - Q(v, q 12 )Q(w, q 2i ) 



Q(v - w : q 12 ) 



which are reduced to ( TSj) after applying first (O) to the left-hand side. 



5 Discussion 

We have calculated the dynamical r matrix for the elliptic Calogero-Moser model 
containing the spectral parameter. We have shown that the Poisson algebra of L and 
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r is not closed: r satisfies a generalized version of the classical Yang-Baxter equation 
(|36p which involves a new object — the X matrix. In fact, there might exist more 
Yang-Baxter type identities involving for instance {X^ 123 \ L^} etc. which might 
give rise to an infinite sequence of r matrices of higher dimensions. A remarkable 
feature of the matrices r and X for the elliptic Calogero-Moser model is that they 
depend only on coordinates q a , cf. 0. Hopefully, this property will hold also for the 
higher r matrices. 

All these questions are still waiting for an investigation. To resolve them, proba- 
bly some geometrical interpretation of the L operator and r matrics might be useful. 
For the case of trigonometric potential and the L operator containing no spectral 
parameter such an interpretation in terms of the hamiltonian reduction exists |j[ [L6| 
and has been applied recently to construct geometrically the corresponding r matrix 



ITf . For the cases of elliptic potential and the L operator with spectral parameter 
no geometrical interpretation is known so far. 

The dynamical r matrix QTSP could provide, possibly, a mean to construct a 
separation of variables for the Calogero-Moser model in the same manner as in the 
case of numerical r matrices for integrable magnetic chains |T7|1 . The work in this 
direction is in progress. 

About the quantum version of the elliptic Calogero-Moser model almost nothing 
is known so far, except a construction of commutative quantum integrals of motion 
A possible way to approach quantization of the model could be to look for a 



quantum version of the relations (||) and 

A natural generalization of the Calogero-Moser model is its relativistic version 



T9| . In this case, the Poisson bracket of L operators can also be represented in an 
r matrix form, though the right-hand-side is now quadratic in L. The r matrices 
again depend only on the coordinates q a , see jI2J and [[2tJ for the cases, respectively, 



of the L operator without, and with the spectral parameter. 
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